Observational constraints on dark energy with a fast varying equation of state 
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We place observational constraints on models with the late-time cosmic acceleration based on a 
number of parametrizations allowing fast transitions for the equation of state of dark energy. In 
addition to the model of Linder and Huterer where the dark energy equation of state w monotonically 
grows or decreases in time, we propose two new parametrizations in which w has an extremum. We 
carry out the likelihood analysis with the three parametrizations by using the observational data 
of supernovae type la, cosmic microwave background, and baryon acoustic oscillations. Although 
the transient cosmic acceleration models with fast transitions can give rise to the total chi square 
smaller than that in the A-Cold-Dark-Matter (ACDM) model, these models are not favored over 
ACDM when one uses the Akaike information criterion which penalizes the extra degrees of freedom 
present in the parametrizations. 



I. INTRODUCTION 

The discovery of the late-time cosmic acceleration Qj 
opened up a new research arena for cosmologists, as- 
trophysicists, and particle physicists. From the view- 
point of particle physics the cosmological constant nat- 
urally appears as a vacuum energy of quantum fields, 
but its energy scale is usually very different from the ob- 
served dark energy scale [2J]. As an alternative to the 
cosmological constant, dynamical dark energy models- 
such as quintessence @, k-essence Q, f(R) gravity J5J], 
f[R,Q) gravity @, DGP braneworld Q|, and Galileon 
[8|-have been proposed. These models give rise to a time- 
varying equation of state w(a) of dark energy, where a 
is the scale factor in the Fricdmann-Lemaitre- Robertson- 
Walker (FLRW) cosmological background |9| . 

At the background level it is possible to discriminate 
between a host of dark energy models by confronting w(a) 
with the observations of Supernova type la (SN la), Cos- 
mic Microwave Background (CMB), and Baryon Acous- 
tic Oscillations (BAO). For this purpose several different 
parametrizations of w(a) have been proposed-which are 
mostly based on two parameters wq and w\ |10l4l5l | (see 
Refs. [la . [l7| for the parametrizations of the Hubble pa- 
rameter H or the luminosity distance Dl instead of w). 
A well known example is the so-called Chevalier-Polarski- 
Linder (CPL) parametrization w(a) = wq + u>i(l — a) 
[HI LUl, where wq is the value of w today (a = 1). The 
2-parameter parametrizations have been widely used for 
constraining the property of dark energy [13 - 1201 ] . 

With two parameters one usually fixes wq and the value 
of w in the asymptotic past (= w p ). In this case it is gen- 
erally difficult to accommodate the time at and the width 
r of the transition in two asymptotic regimes. Bassett et 
al. |2l( first proposed a 4-parameter parametrization in- 
volving at and r. Corasaniti and Copeland [22[ further 



developed this issue and proposed a kink parametriza- 
tion given by w(a) = w + (w p — wq)[(1 + e a */ T )(l — 
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This allows for 



quintessence models with tracker solutions having a rapid 
transition [23|, |24j , which is difficult to be addressed by 
the 2-parameter parametrization. 

Bassett et al. [25[ carried out the likelihood anal- 
ysis with the kink parametrization by using the Gold 
SN la data [1| in 2004. They found that the best-fit 
corresponds to the case in which w is nearly constant 
(w ~ w p — —0.41) for the redshift z larger than 0.1 and 
rapidly decreases toward wq ~ —2.85 for z < 0.1. This 
evolution of w is outside the limits of the two-parameter 
parametrizations, which implies that two parameters are 
not generally sufficient to implement such a rapid transi- 
tion 1 . Corasaniti et al. [331 also showed that the rapidly 
varying equation of state is consistent with the CMB and 
large-scale structure data accumulated by 2004. 

For the kink parametrization the Hubble parameter H 
cannot be derived analytically in terms of a function of 
a. Instead Linder and Huterer (LH) |3J| proposed the 
4-parameter parametrization w(a) = Wf + (w p — Wf)[l + 
(a/at) 1 / r ]~ 1 , which also allows a rapid transition (where 
Wf is the value of w in the asymptotic future). In this 
case there exists an explicit integrated form of H{a) with 
respect to a, so it is technically convenient. Moreover this 
parametrization can accommodate tracker scaling solu- 
tions (w p = 0) with the rapid decrease of w [22|, [23[ and 
thawing quintessence models (w p = —1) with the fast 



In some of quintessence and k-essence models such as thawing 
and tracker models, it is possible to derive the analytic forms 
of w(a) approximately I28H - I321 . Apart from the tracker models 
with the inverse power-law potential 1321 . the field equation of 
state usually contains more than 3 free parameters [2343111 . 



growth of w [35|, |36( . 

For both the kink and the LH parametrizations the 
dark energy equation of state either increases or de- 
creases monotonically. Meanwhile there are some models 
in which w has a minimum-such as quintessence [23|, 
f(R) gravity |23|, and coupled dark energy |38| models. 
In order to implement models in which w has an ex- 
ternum, we propose two new parametrizations given in 
Eqs. ([7]) and (|12l) below. These are based on four param- 
eters at, r, Wp, and wq, which allow fast transitions of w. 
Moreover, in both cases, there exists analytic expression 
of the Hubble parameter 2 . 

In this paper we shall place constraints on the model 
parameters of the LH parametrization ((5|) as well as those 
of the parametrizations ([7]) and (|12[) by using the recent 
observational data of SN la, CMB, and BAO. In each 
model the five parameters at, t, w p , wq (or Wf), and ttm 
(today's density parameter of non-relativistic matter) are 
varied in the likelihood analysis. We also carry out the 
4-parameter space analysis by fixing wo with a number 
of different values between —1 and 0. In order to accom- 
modate the thawing-type models with fast transitions, 
we shall further set w p = — 1 and study the viability of 
the models |Q and (JT2"jl (including transient acceleration 
models) in the 3-parameter space. Note that the observa- 
tional constraints on kink-likc parametrizations different 
from those mentioned above (like those based on the de- 
celeration parameter q) have been studied by a number 
of authors g3H|- 

This paper is organized as follows. In Scc.|H]wc present 
the three parametrizations of w(a) as well as the corre- 
sponding Hubble parameter H(a). In Sec. IIIII we show 
the method of our likelihood analysis to confront the 
models with observations. In Sees. IIV1 IVl IVII we place 
observational constraints on the model parameters of the 
parametrizations (0), (0, and (fT2j) . respectively. Sec. lVIII 
is devoted to conclusions. 



II. PARAMETRIZATIONS OF DARK ENERGY 



written in an integrated form 



/Ode (a) 



(o) 
Pde cx P 



— (1 + w) da 



(2) 



where pj^g is the dark energy density today (a = 1). 
The energy density of non-relativistic matter is given by 

Pm{o) = pm a~ 3 , where p m is its today's value. 
The Friedmann equation gives 



3H 2 = 8TrG{ Pm + p DE ) 



(3) 



where G is the gravitational constant. This can be writ- 
ten as 



H 2 {a) 

Hi 



ngpo-' + ci-ng?) 



exp 



/> 

Ja a 



w)da 



(4) 

where Hq is the present value of H, fl\ n = 
8tt Gpfn /(SHq), and we used the fact that 

Next, we study the parametrization of dark energy al- 
lowing fast evolution of w. One of the examples is given 



, . W p - Wf 
W{a) = Wf + — — y- 

1 1 + (a/a t ) 1/T 



(Model 1) 



(5) 



where a t (> 0) is the scale factor at the transition epoch, 
and t (> 0) characterizes the width of the transition. 
In the asymptotic past (a — ► 0) and future (a — > oo) 
one has w — > w p and w — > Wf, respectively. For the 
parametrization ([5]) the r.h.s. of Eq. (JH) is integrated to 
give 



H\a) 



nW«- 3 + (i-oW) 



/a 



-3(l+u) p ) 



1/t , !/t\ 3T ( W P- W f) 
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l + a t 



l/r 



(6) 



We consider the flat FLRW background described by 
the line element ds 2 = —dt 2 +a 2 (t)dx 2 , where t is cosmic 
time. We take into account dark energy with the time- 
varying equation of state w(a) and non-relativistic mat- 
ter with the density parameter fi„i today. We assume 
that the dark energy density prjE satisfies the continuity 
equation 



Pde + 3H(1 + w)pue = , 



(1) 



where a dot represents a derivative with respect to t, and 
H = a/a is the Hubble parameter. This equation can be 



which is convenient in confronting the model with obser- 
vations. 

However, the parametrization ([5]) does not accommo- 
date the models in which w has an extremum. In order to 
address such cases, we propose the following parametriza- 
tion 



w(a) = w p + (wq - w p ) 



o[l - (a/a t ) lA 



f 



-1/t 



(Model 2) , 

(7) 
where a t > 0, r > 0, and w p ,wq are the values of w 
in the asymptotic past and today, respectively. For the 
parametrization (J7J the Hubble parameter can be ex- 
pressed as 



2 When the Hubble parameter is analytically available, one can 
also determine the Om diagnostic introduced in Ref. |39|| . 



H 2 (a) 
Hi 



= a ( °v 3 



(l-OW)a- 3 ( 1+ ^)exp[/(a)], (8) 



where 



III. DATA ANALYSIS 



/(a) = 3{w - w p ) 
„ ! + (! 



-l/rx 



o{[(o/o t ) 



l/r 



l]r-l} 



(l-oT 1/T )(l+r) 
The equation of state ([7]) has an extremum at 



with the value 



w(a*) = uip + 



(w - w p )t t (t + 1) T l a t 



1 



-l/r 



.(9) 



(10) 



(11) 



If < at < 1 and iu p < Wo, or, a t > 1 and ui p > Ufa, 
then w has a minimum at a = a,. On the other hand, if 
< a t < 1 and w p > wq, or, at > 1 and w p < wo, w has 
a maximum at a = a*. 

For the models characterized by w p > wq with a min- 
imum of w at < a* < 1 (such as quintessence mod- 
els in Rcf. [23|), the transition redshift needs to satisfy 
the condition a t > 1. From Eq. ([IT)]) it follows that 
Ot/e < a* < at and hence a* > 1/e. This means that, for 
w p > wo, the parametrization ([7]) does not accommodate 
the case in which w has a minimum at low redshifts. In 
order to improve this shortcoming, we shall also consider 
the following parametrization 

V /T [l-(a/a t ) 1/T l 
w(a) =Wp + (w -w p ) - M^ — - (Model 3) , 



1 — a t 

where at > and r > 0. Then w has an extremum at 

a t 



(12) 
it 

(13) 



with the value 



w(a*) =w p + - 



1 (iuq - WpK 



1 



-l/r 



(14) 



The equation of state has a minimum either for < at < 
1 and Wp < Wo, or, for a t > 1 and w p > iuo- Since 
a* — > for t>1, we can cover the case of small a* even 
for a t > 1 and w p > u>o- 

The Hubble parameter corresponding to the 
parametrization (TT2"]) is given by Eq. ((SJ), where the 
function /(a) is 



f(a) = 3(w - w P )i" 



,-Vr 



■a 1 /-[(a/a t ) 1 / T -2] 



2(l-a t - 1/T ) 



(15) 



In the regime < a < 1 the parametrizations (0, (|7|), 
and (|12[) can recover the CPL parametrization w(a) = 
wo + w»i(l — a) in the limit that a t 3> 1 (with r = 1 for 
Model 1 and Model 3). 



In this section we explain the method employed to con- 
strain Models 1,2, and 3 observationally. In our analysis 
we use the three datasets: 1) the SN la (Constitution 
H|); 2) the CMB shift parameters (WMAP7) [H|; 3) 
the BAO (SDSS7) [H|. The flat Universe is assumed 
throughout the analysis. 

In SN la observations the apparent magnitude m{z) 
at peak brightness is related with the luminosity dis- 
tance dz(z) = (1 + z) f Q i/ _1 (z)dz through m(z) = 
M + 51og 10 (di(z)/10pc), where z = \ja — 1 is the red- 
shift and M is the absolute magnitude [l[ ■ We define the 
distance modulus 

fi{z) = m{z) - M = 5 \og w [Hod L (z)} + [i , (16) 

where no = 42.38 — 51og 10 /i with h = 
-ffo/IlOOkmscc" 1 Mpc^ 1 ]. The chi square associated 
with SN la observations is given by 



N 



XSN 



la - 2^ 



Mobs(^) ~ M( z ») 



(17) 



where N is the number of the SN la datasct, fi bs(zi) are 
the observed values of the distance modulus, and a^i 
are the errors on the data. We employ the Constitution 
dataset with the total of 397 SN la in order to find the 
minimum of (jTTJ) and the corresponding best-fit param- 
eters. 

The position of the CMB acoustic peaks can be quan- 
tified by the following two parameters [43] 



n = \ n 



>(o) 



dz 



= ^^1, (18) 



/„ H(z)/H ' r s (z 

where z* is the redshift at the decoupling epoch, 



da (z*) = 1Z/[Hoy£lm] is the comoving angular diam- 
eter distance to the last scattering surface, and r s (z*) is 
the sound horizon defined by 



r s (z*) 



dz 



* ^(z)v/3{l + 3^ u V[4^(l + z)]} 



)(o), 



(19) 

Here fl^ and Qiy are the density parameters of baryons 
and photons, respectively. For the redshift z* there exists 
the following fitting formula |48| 

a. = 1048 [1 + 0.00124(^ t V)-°- 738 ] [1 + oi (fi^fe 2 ) 92 ], 

(20) 
where 9l = 0.0783 (n< 0) /i 2 )-°- 238 /[l + 39.5 (Tl^h 2 ) - 763 } 

and g 2 = 0.560/[l + 21.1 (Q^/i 2 ) 1 - 81 ]. The chi square for 
the WMAP7 measurement is 



XcMB 



l cmbCcmb^cmb ■ 



(21) 



where X£ MB = (l a - 302.09, 11-1.725, z*- 1091.3), and 
the inverse covariance matrix is given by [19J 



2.305 29.698 -1.333 
:9.698 6825.27 -113.1! 
-1.333 -113.18 3.414 



c cmb = I 29 - 698 6825.27 



(22) 



The BAO observations constrain the ratio tbao(z) = 
fs{zd)/Dv{z), where r s (zd) is the sound horizon at which 
the baryons are released from the Compton drag of pho- 
tons (denoted as the redshift Zd)- Dy(z) is the effective 

BAO distance defined by D v (z) = [d { a c] \z) z/H(z)] 1 / 3 
[49| . where d a c) (z) = $* H- X {z)dz. For the redshift z d 
there is the following fitting formula [50] 



Zd 



1291 (nl° ] h 2 ) - 251 

1 + 0.659 (fi^/i 2 ) - 828 



( Q )h2\b2-\ 



[1 + hWh 2 ) 1 *], (23) 



where 61 = 0.313 {il^ h 2 )~ 0A19 [l + 0.607 (il^h 2 )"- 674 } 

and b 2 = 0.238(0.$ h 2 ) 0223 . The chi square associated 
with the BAO measurement is given by 



Xbaoi 



X 



Crao^bao , 



BAO^BAO 



(24) 



where Xg AO = (r B Ao(0.2)-0.1905,r B Ao(0.35)-0.1097), 
and the inverse covariance matrix is 1461 



r- 1 

u BAO 



30124 -17227 
-17227 86977 



(25) 



We also use the BAO data from the WigglcZ and 
6dFGS surveys. These data are given in terms of A(z), 
where its theoretical value is 



Ath(z) 



D v {z)dnt ] H 2 



(26) 



and the data are A wigglcZ (z = 0.6) = 0.452 ± 0.018 j51| 
and A MFGS {z = 0.106°) = 0.526 ± 0.028 (H- The chi- 
square is given by 



2 2 

2 _ v^ I A ( z i) ~ A th{zi) 

XBA02 — 2-^i 



1=1 






(27) 



Therefore, the total chi-square from the three datascts 



is 



X - XsNIa + XCMB + XBAOl +XBA02 • ( 28 ) 

The best-fit corresponds to the model parameters for 
which the x 2 is minimized. 



IV. OBSERVATIONAL CONSTRAINTS ON 
MODEL 1 




Figure 1: The dark energy equation of state w versus the scale 
factor a for Model 1 with several different model parameters. 
The line (a) corresponds to the best-fit case of Eq. (|29p de- 
rived by varying the 5 parameters w p ,Wf,at,T, £lm in the 
likelihood analysis. The line (b) shows the best-fit derived 
with the priors w p > and at > 0.5. The lines (cl), (c2), 
(c3) represent the best-fits where the 4 parameters w p ,at,r, 
Qm are varied with the present value of the equation of state 
fixed at wo = —1/3, —0.7, — 1, respectively. 



We first vary the 5 parameters w p ,Wf,a t ,T, and flm 
in the likelihood analysis. The priors on each parameter 
are set to be -10 < w p < 10, -10 < Wf < 10, a t > 0, 



r > 0, and 0.15 < 
parameters are 



(0) 



< 0.4. We find that the best-fit 



w p = 0.141866 



w f = -1.02862 , a t = 0.132023 , 
r = 0.360069 , Q$ = 0.290346 , (29) 



with x 2 — 467.77. In Fig. Q] we plot the evolution of w 
for the best-fit case [line (a)]. Initially there is a period 
where w stays nearly constant (w ~ 0.14), which is fol- 
lowed by the decrease of w around the redshift z larger 
than 10. The dark energy equation of state crosses the 
cosmological constant boundary (w = —1) around z = 1 
and it approaches the asymptotic value Wf = —1.028. 

If w starts to evolve from the value larger than in 
the deep matter era, it is required that the transition 
to the regime w « — 1 occurs in the early cosmological 
epoch (for z larger than 1). In fact, if we carry out the 
likelihood analysis with the priors w p > and at > 0.5, 
the best-fit model parameters are found to be 



We place observational constraints on Model 1 accord- 
ing to the method explained in Sec. IIIII 



w p = 0.0120045 , wf = -1.35856 , a* =0.5, 
t = 0.130893 , Q$ = 0.300473 , (30) 
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w p 


a t 
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n {0) 

,1 £, m 


x 2 





-1.04279 


1.20514 


0.00277414 


0.276338 


470.825 


-1/3 


-1.16253 


1.31208 


0.0502325 


0.280804 


471.530 


-0.5 


-1.17027 


1.30959 


0.0659688 


0.281219 


470.996 


-0.6 


-1.17165 


1.26349 


0.0779737 


0.281295 


470.707 


-0.7 


-1.17077 


1.16698 


0.0907055 


0.281248 


470.456 


-0.8 


-1.16655 


1.18425 


0.116396 


0.280961 


470.276 


-1 


-1.15384 


0.742678 


0.15533 


0.279799 


470.387 



Table I: The best-fit model parameters (4 parameters in total) 
and x 2 f° r Model 1 with several given values of wq . 



with x 2 = 505.983. In this case, the bound on at is 
saturated at a t = 0.5 and the best-fit x 2 is much larger 
than that corresponding to Eq. ([29| . Since r <C 1, the 
transition from the regime w > to the regime w ~ — 1 
occurs quite rapidly. In Fig. [1] we compare the behavior 
of the two best-fits in Eqs. (gH) and (J3TI[) . 

The likelihood analysis of Bassett et al. [25| for the 
kink parametrization, with the SN la and CMB data ac- 
cumulated by 2004, showed that the best fit corresponds 
to the fast transition in low redshifts (z < 0.1). However, 
inclusion of the BAO data as well as the more updated 
SN la and CMB data, seems to point towards a much 
earlier transition from the regime w ~ to the regime 
close to w = — 1. 

In order to study the possibility of the late-time transi- 
tion further, we also study the case in which the value of 
w today (= w ) is fixed. Since Wf = a t [wo{\ + a t /T )- 



the parametrization $5§ can be expressed as 



w(a) 



«',; 



,4/r 



W (l 



-1/r 



W v 



1 + {a/a t yh 



(31) 



For several given values of wq we vary the 4 parameters 
w p , a t ,r, and £l m with the priors —10 < w p < 10, a t > 0, 
t > 0, and 0.15 < fi^ ) < 0.4. In Table Q] the best-fit 
model parameters and the corresponding \ 2 are shown 
for w = 0,-1/3,-0.5,-0.6,-0.7,-0.8,-1. In Fig. Q] 
we also plot w versus a for several different best-fit cases 
{w = -1/3, -0.7,-1). 

For the values of wo between —1 and 0, the initial 
evolution of w for each case shown in Table U exhibits 
a common property. The dark energy equation of state 
is nearly constant with w less than —1 during the deep 
matter era, which is followed by the growth of w in the 
low-redshift regime (z < 1). The parameter r tends to 
be smaller for larger wq, so that the transition becomes 
sharper. This property can be confirmed by comparing 
the three best-fit cases (cl)-(c3) in Fig. Q] 

In Table U we find that \ 2 is more or less similar for 
different choices of w$ between —1 and 0. The best- 
fit ACDM model corresponds to tt$ = 0.269431 with 
X 2 = 471.89, whose x 2 is larger than those given in Ta- 
ble UJ This implies that the transient cosmic acceleration 
models with rapid transitions of w are not excluded by 
the current observational data. 



We need to caution, however, that the parametriza- 
tion ([3"Tj) with given wo has 4 free parameters to fit the 
models with the data, while the ACDM model has only 
one free parameter (ilm ). In order to compare the mod- 
els with different number of free parameters, we employ 
the Akaikc Information Criterion (AIC) [11]. The AIC is 
defined as 



AIC 



Amin 



2V, 



(32) 



where Xmin i s ^ ne minimum value of x 2 1 an d ^ is the 
number of free parameters for each model. For smaller 
AIC the model is more favored. If the difference of AIC 
between two models is in the range < A (AIC) < 2, 
the models arc considered to be equivalent. On the other 
hand, if A(AIC) > 2, one model is favored over another 
one. 

The flat ACDM model corresponds to AIC = 473.89, 
whereas the transient acceleration models in Table U give 
rise to larger values of AIC (e.g., AIC = 478.825 for u> = 
0). The best-fit case ([29)) with 5 parameters corresponds 
to AIC = 477.77. According to the AIC, Model 1 with 
5 or 4 free parameters is not favored over the ACDM 
model. 



V. OBSERVATIONAL CONSTRAINTS ON 
MODEL 2 

Let us proceed to observational constraints on Model 
2. We first vary the 5 parameters w p ,wo,at,T, and fl, n 
in the likelihood analysis. We set the priors on each pa- 
rameter, as — 10 < w p < 10, — 10 < wq < 10, at > 0, 



r > 0, and 0.15 < O 
arc found to be 



(0) 



< 0.4. The best-fit parameters 



w p = -1.10237 



wq = -0.906508 . 
r = 0.505998 , £1$ = 0.280583 , 



«, 



0.739325, 
(33) 



with x 2 = 470.241. In this case x 2 is smaller than that 
in the ACDM model, but it is larger than that in the 
best-fit case ([29]) of Model 1. 

In Fig. [2] the evolution of w for the parameters ([33[) 
is plotted as the solid line (a). As we showed in Sec. [ill 
w has a minimum at a* given in Eq. (JTUJ) either for (i) 
< at < 1, w p < wo, or (ii) a t > 1, w p > wq. The best- 
fit model parameters ([3"3")) correspond to the case (i) with 
a* = 0.43. The equation of state starts from a phantom 
value w p = —1.10237, which is followed by mild decrease 
of w. For a > a* it starts to increase and reaches the 
present value wo = —0.906508. 

The above behavior of w is different from that for the 
best-fit Model 1 with 5 parameters varied. As we see 
in Fig. [1] the best-fit parameters in Model 1 satisfy the 
condition w p > wo, but in this case Model 2 gives rise 
to a minimum only for o, t > 1. Since a* is limited in the 
range a t /e < a* < at and also r is required to be large 
(t > 1) to have small a*, it becomes more difficult to fit 



-0.50 




a 



Figure 2: The dark energy equation of state w versus a for 
Model 2. The line (a) corresponds to the 5-parameter best- 
fit case given in Eq. (|33p . The lines (bl) and (b2) show the 
best-fits where the 4 parameters w p , a t , r, f2m are varied with 
wo = 0, —0.7, respectively. 
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Table II: The best-fit model parameters (4 parameters in to- 
tal) and x 2 f° r Model 2 with several given values of too- 
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w p 
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x 2 





-1 


0.940663 


0.0319267 


0.277496 


472.984 


-1/3 


-1 


0.916762 


0.0689972 


0.280277 


471.702 


-0.5 


-1 


0.909223 


0.0938171 


0.280760 


471.074 


-0.7 


-1 


0.901454 


0.1517700 


0.280844 


470.476 


-0.9 


-1 


0.897160 


0.7207540 


0.280662 


470.245 



Table III: The best-fit model parameters (3 parameters in 
total) and x 2 f° r Model 2 with w p — — 1 and several given 
values of Wq. 
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Figure 3: The la (inside) and 2a (outside) likelihood contours 
in the (at, r) plane derived by varying the 2 parameters at and 
r with Wp = -1, too = -1/3, and Q.<£> = 0.280277 for Model 
2. The black point corresponds to the best-fit case. 



w with the observational data for at > 1 and w p > wq. 
If w p > wq and < a t < 1, w has a maximum at a = a*. 
However, such cases are also difficult to be compatible 
with the observational data. 

By choosing several different values of u>o (= 
0, —1/3, —0.5, —0.7, —0.9), we also vary the 4 parame- 
ters w p , a t , r, f2, rl in the likelihood analysis. The pri- 
ors are set to be — 10 < w p < 10, a t > 0, r > 0, and 
0.15 < £1$ < 0.4. In Table. [0] we show the best-fit pa- 
rameters and x 2 f° r each wo- In all cases we find that 
< a t < 1 and w p < wq, so that w has a minimum at 
< a* < 1. 

In Fig. [5] we show the variation of w for the best-fit 
cases with wq = and wq = —0.7 as the lines (bl) and 
(b2), respectively. The growth of w in the regime a > a* 
is sharper for larger values of wo. This reflects the fact 
that, in Table [TTl t gets smaller for wq increased. We also 
note that the models with larger wq tend to be disfavored 



because of the increase of x 2 seen i n Table [TTJ 

In addition to wq, we also fix w p to be —1 and vary 
the 3 parameters a t ,r, flm . In Table Hill we present the 
best-fit model parameters for several different choices of 
wo (> —1). In all cases the transition scale factor is in 
the range < a t < 1, so that w has a minimum for 
< a* < 1. In fact, the evolution of w starting from 
— 1 and having a minimum by today is present for dark 
energy models based on f(R) theories [37| (although if 
does not continuously grow for a > 1). Table Hill shows 
that, for larger wq, t tends to be smaller, whereas x 2 gets 
larger. 

In Fig. [3] we illustrate the la and 2a observational con- 
tours in the (a t ,r) plane for w p = —1, Wq = —1/3, and 

flm = 0.280277. This is the marginal case in which the 
Universe enters the phase of cosmic deceleration today. 
The redshift and the width of the transition are con- 
strained to be 0.87 < a t < 0.99 and < r < 0.18 (68% 
CL). Unless the rapid transition occurs at the redshift 



close to today, the model is not compatible with the ob- 
servational data. For larger r the values of w at a = a* 
start to deviate from —1, so that those cases are more 
difficult to satisfy observational constraints. 

As mentioned in Sec. [TV] the AIC for the flat ACDM 
model is AIC = 473.89. For the 5-parameter best-fit 
case in Eq. (|33| and for the 4-parameter and 3-parameter 
best- fit cases given in Tables [TT| and IIII1 the AIC in each 
model is larger than that in the ACDM model with the 
difference more than 2. Hence the AIC criterion shows 
that the ACDM is generally favored over Model 2. 



VI. OBSERVATIONAL CONSTRAINTS ON 
MODEL 3 

Finally we proceed to observational constraints on 
Model 3. When the 5 parameters w p , wo,at,r 7 ilm are 
varied in the likelihood analysis, we set the same priors 
as those given in Model 2. The best-fit parameters are 
found to be 



w p = -1.10733 



w = -0.897454 . 
t = 0.652107 , Q$ = 0.28053 , 



a i 



0.737871 , 

(34) 



with x 2 = 470.235. 

As we see in Fig. 01 the evolution of w corresponding 
to Eq. (|34p is similar to that for the best-fit parameters 
(|33jl of Model 2. Since < at < 1 and w p < wq for the 
model parameters (|34|) . w has a minimum at a* = 0.47. 
The difference between Models 2 and 3 is that even for 
a t > 1 and w p > wq the equation of state for Model 3 can 
take minima with smaller values of a* given in Eq. (1131) . 
However we find that the models with at > 1 and t>1 
are disfavored because w(a*) deviates from — 1. 

For the 4-parameter parametrization with a number 
of different values of wo (= 0, —1/3, —0.5, —0.7, —0.9) we 
also vary the parameters w p ,at, t, ilm with the same pri- 
ors used for Model 2. In Table IIVI we summarize the 
best- fit parameters as well as the \ 2 f° r each wq. For 
wo = — 1/3, —0.5, —0.7 one has < a t < 1 and w p < wq, 
in which cases w has minima at < a* < 1. If wo > —0.5, 
the growth of w in the regime a > a* is very rapid (see 
the line (bl) in Fig. gj). 

For the best-fit parameters corresponding to wq = 
0,-0.9 one has at > 1 and w p < wq. In those cases w 
has maxima at a* larger than 1 and hence w is a growing 
function for a < 1. Since r is extremely small for wq = 0, 
the transition of w occurs almost like a step function. 
However such an instant transition cannot be regarded 
as a realistic model of dark energy. For wq = —0.9, w 
has a maximum (w(a*) = 5.5) at a* = 2.4. In this case 
the evolution of w is not very different from that for the 
best-fit case (|55|). apart from the fact that for wo = —0.9 
the equation of state is a growing function in the regime 
a < 1. 

We also vary the 3 parameters at, t, fim by fixing w p to 




Figure 4: The dark energy equation of state w versus a for 
Model 3. The line (a) represents the 5-parameter best-fit case 
given in Eq. (|34|1 . whereas the lines (bl) and (b2) correspond 
to the best-fits derived by varying the 4 parameters w p ,a,t,T, 
Q.m with u>o = —1/3, —0.9, respectively. 
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0.22092 
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Table IV: The best-fit model parameters (4 parameters in 
total) and \ 2 f° r Model 3 with several given values of wq. 
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0.277967 
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472.013 


-0.5 
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0.904406 
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-0.7 
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0.280590 
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-0.9 


-1 


0.901692 


0.788501 
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Table V: The best-fit model parameters (3 parameters in to- 
tal) and x f° r Model 3 with w p — — 1 and several given values 
of WQ. 



IVlwe show the best-fit values as well as x 2 for each wq. In 
most cases the transition redshifts are around at = 0.9. 
As we increase wo the parameter r gets smaller, so that 
the transition occurs more rapidly. For larger wq, x 2 
tends to be larger. 

In Fig. [5] we plot observational bounds in the (a t ,T) 



be —1 for several different choices of wq (> —1). In Table plane for w p = — 1, wq = —1/3, and Q 
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Figure 5: The la (inside) and 2cr (outside) likelihood contours 
in the (at, r) plane derived by varying the 2 parameters at and 
r with w p = -1, wo = -1/3, and 0$ = 0.279820 for Model 
3. The black point corresponds to the best-fit case. 



Comparing it to Fig. [31 we find that the larger values 
of t can be allowed in Model 3. This reflects the fact 
that in Model 3 the values of w(a*) do not deviate from 
— 1 significantly for r < 0.5. The two parameters are 
constrained to be 0.87 < a t < 0.95 and 0.12 < r < 0.44 
(68% CL). 

For all the best-fit cases discussed above, the AIC is 
larger than that in the flat ACDM model. Hence Model 
3 is not favored over the ACDM model according to the 
AIC. 



VII. CONCLUSIONS 

In this paper we placed observational constraints on 
the three models allowing fast transitions of w, by using 
the data of SN la, CMB shift parameters, and BAO. 
Unlike the 2-parameter paramctrization such as w(a) — 
Wo + Wi(l — a) , the parametrizations ©, ([7]), and (fT2"]l 
have two more parameters at and r by which the time 
and the width of the transition can be accommodated. In 
Model 1 the dark energy equation of state monotonically 
increases or decreases in time, whereas in Models 2 and 
3 w has either a minimum or a maximum depending on 
the values of w p , wq, and at- For all these models the 
Hubble parameter H is analytically known in terms of 
functions of a. 

When the 5 parameters w p ,Wf,a t , r, Qm are varied in 
Model 1, the best- fit parameters are given by Eq. (f2"9"]l 
with x 2 = 467.77. This corresponds to the solid curve 
(a) in Fig. [TJ in which case the equation of state enters 



the regime w ~ —1 in the early cosmological epoch. If we 
put the prior on the transition redshift as at > 0.5, x 2 be- 
comes significantly larger than that without a prior on a t . 
This means that the late-time transition (at > 0.5) from 
the regime w ~ to the regime w ~ — 1 is disfavored ob- 
servationally. If we vary the 4 parameters Wp,at,T,Sl m 
with several different values of wq between —1 and 0, the 
parameter r tends to be smaller for increasing wq- Al- 
though the x 2 hi Model 1 with 5 or 4 parameters can be 
smaller than that in the ACDM model, the AIC shows 
that Model 1 is not favored over the ACDM model. 

The best-fit parameters for Model 2 corresponds to the 
case in which w starts from a phantom value w p = —1.10, 
takes a minimum —1.17 at a* = 0.43, and grows to the 
value wq = —0.91 by today. This is different from the 
evolution of w for the best-fit parameters of Model 1. 
This difference mainly comes from the fact that, in the 
cases w p > wq for Model 2, w has a minimum at a* 
given by Eq. (fTO]) only for a t > 1. While Model 2 can 
accommodate the late-time transition having a minimum 
of w, it is difficult to address the early sharp transition 
with w p > wq. The 4-parameter likelihood analysis for a 
number of fixed wo (between — 1 and 0) leads to similar 
best-fit evolution of w to that for the 5-parameter best-fit 
case, with a faster transition for larger wo. 

In Model 3 the equation of state has an extremum at 
a* = a t /2 T , which can be smaller that that for Model 
2. If w p > wq, however, the early transition of w with a 
minimum requires the condition t>1. This value of r is 
too large to accommodate the early transition compatible 
with observations, because the minimum value of w tends 
to deviate from — 1. The 5-parameter likelihood analysis 
shows that the best-fit case in Model 3 is similar to that 
in Model 2. The likelihood results for 4 parameters (wq 
fixed) and for 3 parameters (wq and w p fixed) also give 
rise to the similar results to those found in Model 2. 

In Models 2 and 3 the AIC is always larger than that 
in the ACDM model with the difference more than 2. 
Hence the models with the late-time fast transition to 
the non-accelerating Universe are disfavored compared 
to the ACDM model. The joint data analysis based on 
SN la, CMB, and BAO prefers the models in which w 
do not deviate significantly from —1 in the low-redshift 
regime. 

Although the minimum value of % 2 in Model 1 with 
5 parameters is smaller than those in Models 2 and 3, 
the minima in Models 2 and 3 are still inside the la 
region corresponding to Model 1. Therefore we cannot 
prefer/exclude any paramctrization with respect to any 
other one. 

Recently it was shown that, in the framework of the 
CPL paramctrization, the observational constraints on 
dark energy are sensitive to the presence of the cos- 
mic curvature fl K [5J]. They found that the CPL 
paramctrization is not sufficiently flexible to model the 
rapidly varying equation of state in low redshifts even 
for fl K =/= (see also Ref . [43[ ) . It will be of interest to 
study how the effect of the cosmic curvature affects the 



observational constraints on the models discussed in this 
paper. We leave this for future work. 
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